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Area  or  constant  defined  in  Eq.  28 
Sonic  velocity 
Constant  defined  in  Eq.  29 


2 2 
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Plateau  pressure  coefficient 
Separation-point  pressure  coefficient 
Specific  heat  at  constant  pressure 
Error  function  as  defined  in  Eq.  37 
Force 
Enthalpy 

Mass  integral  defined  in  Eq.  47 
Momentum  integral  defined  in  Eq.  48 

Separation-point  distance  measured  from  the  nozzle  throat 
Mach  number 

Undisturbed  Mach  number  before  boundary-layer  separation 
Mass  flow  rate 

Velocity  profile  parameter  (see  Eq.  2) 

Pressure 

Heat  transfer  rate  into  the  wall  in  the  separated-flow  region 
Gas  constant 

Reynolds  number  of  the  undisturbed  flow  before  the  separation 
point 


» 

h 


NWC  TP  5788 


Temperature 

Velocities  in  x,  y direction,  respectively 
Total  velocity 

Coordinates  of  the  reference  system 
Coordinates  of  the  intrinsic  system 
Nozzle  half-angle 
Effective  TVC  deflection  angle 

Dummy  variable  associated  with  the  error  function 
Specific  heat  ratio 
Boundary  layer  thickness 

Height  of  the  subsonic  portion  of  the  boundary  layer 
Energy  similarity  parameter 
Nozzle  expansion  ratio,  A/A* 

Dimensionless  coordinate,  <r y/x 
Momentum  similarity  parameter 
Angle 

Prandtl-Meyer  expansion  angle 
Shock  angle 

Slip-line  angle  (Fig.  21) 

Total  enthalpy  ratio  of  the  jet,  h /h 

J o oe 

Molar  specific  heat  ratio  of  the  jet  defined  in  Eq.  39  or 
\ i /2  as  defined  following  Eq.  7 

Mass  similarity  parameter 

Molecular  weight 

Prandtl-Meyer  angle 

Dimensionless  quantity,  y/6 
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Density 

Jet  spreading  parameter 
Shear  stress 

Velocity  ratio,  u/u  or  u/u 
c e 

Temperature  ratio,  T /T 
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Shear  work  and  heat  transfer  across  the  mixing  layer 
Energy  transfer  by  the  entrained  mass  flow 
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Sonic  conditions 

Attached  side 
Separated-flow  region 
Separated  point 
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Shear  layer  edge 
Gas 

Separating  streamline  or  TVC  gas 

Coordinate  shift  in  the  mixing  layer  theory  due  to 
momentum  conservation 

Stagnation  or  undisturbed  flow 

Port  side  or  particle 

Subsonic  portion  of  the  boundary  layer 

Wall 

Components  in  x and  y directions,  respectively 
Free  stream 
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INTRODUCTION 


The  Navy's  off-axis  attack  missile  concept,  currently  being 
developed  by  the  Naval  Weapons  Center  ( NWC),  China  Lake,  is  directed 
toward  a low-cost,  high-density,  maneuverable,  air-to-ground  weapon 
for  attack  aircraft  with  single-pass  capability.  Design-to- cost  consi- 
derations dictate  use  of  an  inexpensive  control  system  that  is  nonethe- 
less capable  of  providing  rapid  response  and  high-incidence  airframe 
maneuvers.  These  considerations  have  caused  the  Navy  to  investigate 
a boundary-layer  thrust  vector  control  (BLTVC)  system. 

The  fundamental  principle  of  BLTVC  operation  involves  local 
modulation  of  the  so-called  "back  pressure"  about  the  flow-separation 
region  in  a highly  overexpanded  rocket  exhaust  nozzle.  This  modula- 
tion is  done  by  selectively  allowing  inflow  of  ambient  air  into  the  back- 
pressure region,  which  in  the  unperturbed  state  is  at  low  pressure 
because  of  viscous  pumping.  Natural,  unforced  flow  of  ambient- 
pressure  air  through  control  ports  into  the  back-pressure  region  causes 
the  separated,  supersonic  rocket  exhaust  flow  to  reattach  to  the  nozzle 
on  the  side  opposite  controlled  inflow.  Since  the  reattached  flow  thence 
follows  only  a portion  of  the  nozzle  wall,  the  rocket  exhaust  can  be 
vectored  at  an  angle  that  approaches  the  nozzle  half-angle,  thus  pro- 
viding a side  force  consistent  with  the  lateral  exhaust  momentum  (also 
manifested  as  an  asymmetrical  nozzle  pressure  distribution).  Recent 
Navy-funded  work  by  Chandler  Evans,  Inc.  (Ref.  1 and  2)  has  resulted 
in  development  of  a viable  BLTVC  hardware  concept.  The  phenom- 
enology of  this  concept  has  been  variously  called  fluidics,  boundary- 
layer  control,  and  the  Coanda  effect  since  about  1930.  References  3, 

4,  and  5 are  indicative  of  the  historical  background  available. 

Although  analytical  studies  in  related  technology  areas  have  been 
conducted  by  McDonnell  Douglas  Astronautics  Company  (MDAC)  and 
others,  no  theoretical  studies  before  the  present  one  have  directly 
addressed  the  evaluation  of  BLTVC  performance.  The  successful  con- 
duct of  such  a study  offered  the  twofold  advantage  of  providing  infor- 
mation for  predicting  system  performance  and  also  of  providing  a 
basis  for  critical  evaluation  of  experimental  programs.  With  these 
considerations  in  mind,  the  present  study  encompassed  the  following 
objectives:  (1)  development  of  analytical  or  semiempirical  descrip- 
tions of  the  dominant  phneomena  that  characterize  the  thrust  vector  con- 
trol (TVC)  flow  field,  (2)  derivation  of  similarity  parameters  that  affect 
TVC  nozzle  performance,  (3)  development  of  preliminary  guidelines 
for  flight-environment  simulation  with  ground  experiments,  and 
(4)  generation  of  design  charts  for  use  in  TVC  performance  evaluation. 
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This  report  discusses  overall  flow-field  conditions,  presents  most 
of  the  analyses  of  flow-field  modeling,  and  gives  computational  results. 
BLTVC  similarity  and  scaling  studies  are  described;  guidelines  are 
given  for  nozzle-design  matters  such  as  port  location  and  size,  nozzle 
expansion  ratio,  and  chamber -pressure  range;  and  conclusions  and 
recommendations  are  summarized. 

It  should  be  noted  that  this  study  used  results  of  cold-flow  and 
static  rocket  tests  exclusively.  These  tests  were  conducted  in  a 
quiescent  environment,  and  caution  should  therefore  be  exercised 
regarding  the  applicability  of  these  data  for  direct  evaluation  of  flight 
performance  without  consideration  of  the  flow  phenomena  that  are 
peculiar  to  BLTVC  operating  environments. 
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FLOW-FIELD  CONDITIONS 

To  date,  all  BLTVC  cold-flow  and  rocket  tests  have  been  conducted 
without  ambient  flow  external  to  the  nozzle.  The  difference  in  BLTVC 
performance  between  a cold-flow  or  rocket  static  test  and  a wind- 
tunnel  or  flight  test  may  be  very  significant  because  of  the  presence  of 
the  ambient  flow.  As  shown  schematically  in  Fig.  1,  the  main  differ- 
ence between  a BLTVC  nozzle  with  and  without  ambient  flow  is  the 
interacting  Bow  field,  which  involves  a two-stream  jet-mixing  process. 


BLTVC  IN  A QUIESCENT  EXTERNAL  ENVIRONMENT 


For  BLTVC  nozzles  without  external  ambient  flow,  the  injection 
mass  does  not  interact  strongly  with  the  nozzle  main  flow  in  a conven- 
tional sense.  However,  it  does  trigger  the  highly  unstable  separated 
main  flow  of  an  overexpanded  nozzle  to  attach  on  the  side  opposite 
the  injection  ports.  This  main-flow  attachment  process  appears  to  do 
two  things:  ( 1)  it  provides  an  ejector-type  action  that  carries  the 

injection  mass  out  of  the  nozzle  adjacent  to  the  affected  boundary,  and 
(2)  it  creates  a large  volume  in  the  aft  ends  of  the  port  sides.  Conse- 
quently, a significant  feedback  of  the  ambient  air  into  the  inside  rear 
surface  on  the  nozzle  port  side  occurs,  and  the  port-side  pressure 
approaches  the  ambient  pressure  level,  p^ , very  rapidly. 


The  cylindrical  section  of  the  nozzle  shown  in  lig.  1 assists  in  the 
attachment  process.  The  pressure  is  lower  on  the  attached  side  when 
the  flow  does  reattach.  Without  reattachment,  the  pressure,  influ- 
enced by  the  ambient-flow  feedback,  rises  more  rapidly,  and  this 
reduces  the  attainable  side-force  level.  Furthermore,  hysteresis  may 
occur  in  the  side  force  during  the  cycle  of  rising  and  falling  chamber 
pressure,  as  shown  in  Fig.  2.  This  happens  because  flow  attachment 
does  not  occur  at  the  same  p^  value  when  chamber  pressure  increases 
and  decreases.  Recent  cold-flow  tests  employed  a longer,  cylindrical 
aft  section  that  eliminated  the  hysteresis  loop  by  minimizing  the 
ambient-flow  feedback  effect.  Typical  results  are  shown  later  in 
Fig.  34. 


BLTVC  WITH  EXTERNAL  AMBIENT  FLOW 

For  TVC  nozzles  that  operate  with  external  ambient  flow,  the 
mechanism  that  triggers  the  main-flow  attachment  process  is  strongly 
influenced  by  the  free  shear  layer  formed  by  the  ambient  flow.  This 
free  shear  layer  interacts  with  the  free  shear  layer  along  the  boundary 
of  the  main  flow  and  causes  a two-stream  mixing  process,  as  shown 


3 


* 


FIG.  1.  Flow  Field  Conditions  for  BLTVC  Nozzle  With  and  Without  External  Ambient  Flow. 
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FIG.  2.  Hysteresis  Loop  Observed  in  Cold-Flow  Tests. 
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in  Fig.  1.  The  injection  mass  affects  the  flow-field  conditions  in  the 
regxon  bounded  by  the  two  free  shear  layers  and  the  inside  surface  of 
the  port  side.  Therefore,  there  is  no  direct  ambient-flow  feedback, 
and  the  pressures  on  the  port  side  are  more  uniform  throughout  the’ 
entire  region.  Since  no  tests  that  include  ambient-flow  simulation  have 
been  conducted  to  date,  the  effect  of  ambient  flow  on  overall  BLTVC 

Pf  ^or^ance  must  be  estimated  by  using  a rational  engineering  model 
of  the  flow  interactions. 

• u r°?  a confiSurational  standpoint,  the  cylindrical  aft  section 
wmch  helps  to  achieve  high  side  force  in  cold-flow  and  rocket  static 
tests,  may  not  be  required  for  wind-tunnel  or  flight  tests  because  jet 
mixing  with  ambient  flow  occurs  on  both  the  attached  and  unattached 
sides.  This  process  tends  to  assist  flow  reattachment  and  should  be 
quantitatively  investigated  by  experiment. 


FLOW -REATTACHMENT  PROCESS 


On  the  attached  side,  the  nozzle  main  flow  attaches  to  the  cylindri- 
cal section  through  a recompression  process  similar  to  the  reattach- 
ment process  that  is  characteristic  of  a ramp-induced  separated  flow. 
As  a first  approach,  it  seems  suitable  to  use  the  results  obtained  in 
some  ramp-induced  separation  and  reattachment  cases  to  approximate 
the  flow-field  conditions  on  the  attached  side  of  a TVC  nozzle.  Refer- 
ence 6 provides  such  results,  as  shown  in  Fig.  3,  4,  and  5.  Approxi- 
mate equations  for  the  curves  shown  in  these  figures  are: 

For  laminar  flow, 


C R = 0.  20993  + 

ps  xo  77J  T 


CPPRxo'/4=0- 37098 + 


For  turbulent  flow, 

C R 1/10  = 0.  10912  + 
ps  xo  7 T 


CPPRJ 110  + 


0. 55034 

0.  06027 

(M  -1) 

7 

<M  -ir 

o 

1.  03886 

0. 14985 

(M  -1)  " 

7 

(Mq-Ij 

1.  18214 

0.  21857 

(Mo-l) 

<Mo-l|2 

1.  33474 

0. 35426 

(Mo-l) 

(Mo-l)2 

where  Cps  and  Cpp  are  the  separation-point  and  plateau-pressure 

l°e“1C1tnt8'  re3PectlvelY-  Rxo  is  the  Reynolds  number  and  M„  is 
the  Mach  number.  Both  are  evaluated  at  the  beginning  point  of 
interaction. 
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For  the  purpose  of  expediting  hand  calculations,  two  more 
auxiliary  charts  for  unit-length  Reynolds  number  for  various  Mach 
numbers  are  given  in  Fig.  6 and  7. 

Figures  6 through  12  are  applicable  for  attached-side  evaluation 
for  all  conditions.  Different  results  may  be  obtained,  depending  on 
boundary-layer  transition  phenomena.  To  produce  these  curves,  a 
transition  Reynolds  number  of  10^  has  been  used. 


/ PcATL 


-:  i.1  >.t.,.'l'^4:i.,V  r*J 

% 


ili^!'r®:  ;-^  ■ -$^(!.  ;r4;V;,£& S;‘'t- 


4^ 
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FIG.  10.  Plateau  Presaures  on  the  Attached  Side  for  7 = 1.2  and  Mc  = 3.5,4.0,  and  4.5. 
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FIG.  1 1 . Plateau  Pressures  on  the  Attached  Side  for 


FIG.  1 2.  Deflection  Angles  of  the  Attaching  Flow 
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FLOW-FIELD  MODELS  AND  RESULTS 


SEPARATION  POINT 

Flow-field  separation  is  characteristically  present  in 
overexpanded  nozzles.  Methods  to  predict  the  separation  point  insi 
the  nozzie  inchide  a correlation  technique  developed  by  Kalt  and 
Badal  (Ref  7)  and  an  analytical  method  developed  by  Thompson 
(Kef.  8).  Thompson's  analytical  predictions  do  not  agree  well  with 

s^rlnrs^Kalt  ant^dai  /or 

nne  13)<  Thl®  18  Probably  due  to  some  of  the  basic  assump- 

tions that  were  used  in  his  theoretical  development.  For  example 

tum^rr  n , K [lat'plaiie  flow  model  assumed  that  the  momer 
turn  of  the  fluid  between  the  sonic  line  and  the  wall  boundary  is  redr 

to  zero  by  the  action  of  the  pressure  feedback. 

Thompson's  flow  model 
by  Kampe  (Ref.  9). 
applied  to  the  control  volum 
predictions  of  the  r~ 
equating  the  momentum  flux  of  the  sub 
ated  with  the  immediate 


, which  was  used  in  early  BLTVC  si 
, is  based  on  the  law  of  momentum  conservati 
e shown  in  Fig.  14.  With  this  mode 
separation-point  location  have  been  obtained  b 

isonic  boundary  to  that  ass 

pressure  gradient,  i.  e.  , 


uc  0{y/6),  boundary-layer  flow  velocity  at  distance  y from 
the  wall  7 

6s  - height  of  the  subsonic  portion  of  the  boundary  layer 
p = density  7 

pc  - local  boundary-layer  and  mainstream  pressure  at  the 
separation  point 

Pb  ~ pressure  downstream  of  the  separation  point 

It  sWd  be  noted  that  pb  is  the  vortex  (plateau)  pressure  for  the 
attached-  side  separation  and  is  equal  to  the  base  or  ambient  pres- 
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(Note  that  ug/a  - 1 is  being  used.  ) 
Equation  3 becomes 
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n 2 (1  y t 

2 + n uc  (mJ  = (Pco  -Pc) 


(4) 


But  since  Pc  u< 


^PcMc  , it  is  found  that 
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2 + n 


Ic_  _ /£«,  pA 
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or 
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\ Pqc 
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y - l 
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- 1 > = M‘ 
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Results  from  Eq.  6 are  shown  in  Fig.  13. 


For  axisymmetric  flow  with  temperature  variation,  a more  general 
velocity  profile  is  used  to  replace  the  power-law  profile.  The  velocity- 
profile  family  described  in  Ref.  10  by  Kutateladze  and  Leont'ev  is 
suitable  for  the  analysis. 


u 


*- 
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Tw/Te,  ^ - y/6,  and  n = pressure-gradient 


where  u = velocity.  ^ 

parameter.  Subscripts  wW'e  refer  to  wall  and'noTzTe  In*i7c id" 
flow  conditions,  respectively.  Letting 


h 


k f 1 = (T  /T  )^2 

w e 


...  . ‘ .1  • * r^r 


©yapi 


The  following  inviscid-flow  expressions  are  also  used 


WB 
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where 

X = ( T /T  ) 1/2 
w c 

Other  equations  are: 
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separated-f Low  studies  has  shown  that  this  is  a reasonable  approxi- 
mation of  the  actual  profile.  Both  curves  were  obtained  by  integrating 
the  momentum  flux  within  the  shear  layer  from  the  wall  up  to  the 
point  where  the  local  Mach  number  equals  1.3  5.  Comparisons  with 
several  data  correlations  and  theoretical  predictions,  including 
Thompson's  results  (which  were  obtained  by  terminating  the  integra- 
tion at  a local  Mach  number  of  unity),  are  aLso  shown  in  Fig,  13. 

Based  on  these  two  curves,  local  Mach  number  at  the  separation 
point  can  be  computed  for  various  chamber  pressures.  The  results, 
shown  in  Fig.  15,  are  labeled  Attached  Side  for  conditions  both  with 
and  without  attachment  and  Port  Side  (With  TVC)  for  cases  with  TVC. 
Similar  results  for  cases  with  V = 1.  2 are  shown  in  Fig.  16. 

When  the  main  flow  from  the  nozzle  separates,  it  may  or  may  not 
reattach  again  on  the  side  opposite  to  the  port  side.  The  reattachment 
phenomenon  is  influenced  by  the  magnitude  of  the  TVC  gas  injection. 
Since  the  cylindrical  aft  section  has  a deflection  angle  of  15  degrees 
with  respect  to  the  nozzle  wall,  the  attaching  flow  undergoes  a recom- 
pression process.  The  minimum  pressure  attainable  at  the  separation 
point  can  be  determined  rather  easily  by  a simple  graphical  method. 
The  results,  as  shown  in  Fig.  17,  were  obtained  as  follows: 

1.  Local  Mach  number  and  the  corresponding  local  pressure 
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were  calculated  for  various  chamb 
fied  nozzle  geometry.  For  the  pre 
15-degree  turning  angle  is  used,  a 
dashed  lines  in  Fig.  17. 

The  downstream  pressure  can  be  o 
shock  relations  for  Mc  going  throuj 
angle.  The  pressures  downstream 
shown  as  solid  lines  in  Fig.  17. 

When  p2  is  equal  to  p ^ (14.  7 psia  , 
upstream  pressure,  pc,  represent! 
that  the  flow  can  attain  without  losi 
shown  in  Fig.  17,  and  the  correspc 
shown  in  Fig.  15. 


Particle  - Effect  Considerations 


L,pb-Pc>  MgU(Pb-Pc)  6sJp5BL(pb- 

That  is;  the  effect  of  particles  on  momentum  balance 
mathematically  by  two  parameters,  A and  B.  These 
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* 5 PERCENT 


'"CHOKED 


CHOKED 

PLOW 


0.32  (WITHOUT  TVCt 


MASS  FLOW  RATE,  im \G~\J /ApP 

FIC  ” WBh'TVr  ?“rin* ,h'  Ti* « 

J2S2SST  or  0,"txp,n*d 

Since  choked  flow  occut q a f _ ^ 

uuw  ^^urs  at  a pressure  ratio  nf  r 

cuTse7ear\TeVXTVCat  ^ beginninS  of  TVC  mai 

stream  such  that p T ^ sMft  th< 
j Pc/Poo  — 0.67.  Accordingly  th 

drop  approximately  5%,  as  compared  with  the  ch 
Therefore,  the  choked-flow  rh  value  can  be  used 

poses,  values  of  m for  the  cold-flow  tests  were 
are  given  m Fig.  20.  were 

TWO-STREAM  JET-MIXING  INTERACTION 

The  two- stream  jet-mixing  theory  described 
and  application  of  the  original  jet-mixing  theories 
Chapman  (Ref.  12)  and  Korst  (Ref.  13).  It  is  con 
Chapman- Korst  base-flow  model.  Basic  and  rela 
be  found  in  Ref.  12,  13,  and  14. 

The  Bow  model  consists  of  two  elements:  an 
as  shown  in  Fig.  21a  and  a two- stream  mixing- fl< 
in  F lg.  21b.  6 


Inviscid  Jet  Boundaries 

The  inviscid- edge  conditions  are  needed 
layer  properties.  To  compute  these  conditio 
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FIG.  20.  TVC  Mass  Flow  Rate  for  Three-Port  CEN-1 7 Nozzles  Under  Various  Flow  Conditions. 
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sion 


Wedge  theory  is  used  in  order  to  specify  the  recompressed  inviscid 
flow  properties.  The  pressure  rise  is  determined  by  the  position  of 
the  slipline,  the  oncoming  Mach  number,  and  the  shock-wave  angle  of 
each  stream.  For  example,  for  a selected  value  of  the  slipline  angle, 

9 s,  the  recompression  pressure,.  P4,  is  computed  using  oblique  shock 
relations,  viz. 
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Mixing-Layer  Properties 
To  describe  the  mixi 


o y/x  and  a is  the  jet- spreading  parameter  which,  from 
3 17,  is  given  by 

= 12  (1  + 0.  23  M ) 


erf  rj,  is  defined  as 


Although  the  present  velocity  profile  may  not  exactly  follow  the 
error-function  shape,  it  is  widely  known  that  an  integral  method, 

such  as  the  one  used  here,  is  not  very  sensitive  to  small  changes  ii 
profile  shape. 

The  total  enthalpy  profile  is  obtained  using  the  Crocco  relation 
A = <£  + A (1  - = A.  + M_A  ^ A / 


o'^oe  the  total  enthalpy  ratio  of  the  jet 
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The  density  profile  for  a constant-pressure  mixing  region  is 
given  as 


^<fT 


h -u  /2  c ^ 
oe  e p 

h -u  /Z  pe  e 

o r 


A - <f,2C 


where  \=c  ±A(  / c is  the  molar  specific-heat  ratio  of  the  iet. 

p pe  e r J 

The  thermodynamic  properties  such  as  specific  heat,  Cp,  and 
specific-heat  ratio,  y , are  assumed  to  be  independent  of  temperature 
for  the  purpose  of  simplicity  and  computational  efficiency. 

Important  Streamlines 

There  are  two  important  streamlines  in  a mixing-layer  analysis. 

One  is  called  the  separating  (or  dividing)  streamline,  the  other  the 
reattaching  (or  discriminating)  streamline.  The  latter  is  either 
identical  to  or  related  to  the  former,  depending  on  whether  there  is 
an  injection  or  not.  Therefore,  it  is  necessary  to  locate  the  separating 
streamline  first.  To  do  this,  the  coordinate  shift,  rjm,  should  be 
determined.  Since  the  initial  boundary- layer  effects  have  been 
neglected,  the  "restricted"  mixing  theory  of  Korst  for  a negligible 
initial  boundary  layer  is  directly  applicable  (Ref.  14).  The  X-momentum 
entrained  by  the  recirculating  flow  is  usually  ignored.  Hence, 
X-momentum  conservation  gives  the  following  relation  per  unit  width 
of  the  mixing  layer: 


(pu2Y) 


2 , 

pu  dy 


The  ye  is  large  enough  that  <be  — • 1 . 

The  coordinate  system  shown  in  Fig.  22  indicates  the  following 
relation: 


X as  x 


Y s y " ym  (x)  (42 

where  X-Y  and  x-y  are  the  inviscid  and  intrinsic  coordinate  systems, 
respectively.  Using  the  similarity  coordinates,  y = (x/<x)rj,  Eq.  42 
becomes 

Y = V <”  - "J 

Equation  40  then  becomes 
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FIG.  22.  Coordinate  Systems  and  Nomenclature 


it  results  in 


Using  tlie  similarity  coordinates  and  Eq.  39  and  43,  an  equation  results 
that  allows  the  separating  streamline  to  be  determined. 
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(48) 


where  A and  A are  defined  in  Eq.  38  and  39.  Equation  46  becomes 
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dr,  (49) 


Reattaching  (Discriminating)  Streamline.  The  reattaching  stream- 
line, d,  is  the  other  important  streamline  considered.  This  is  the 
streamline  that  divides  the  mass  flow  in  the  mixing  layer  from  the 
flow,  which  is  turned  back  into  the  base  region  to  be  recirculated 
and/or  vented.  When  there  is  no  injection,  the  d and  j streamlines 
are  identical.  With  injection  present,  the  j and  d streamlines  are 
related  by  the  amount  of  mass  pumping  in  or  out  of  the  recirculation 
region  according  to  the  following  relation: 


I 


V. 


m = ± I pu  dy 

'"j 


(50) 


where  the  negative  sign  refers  to  yd  < yj,  i.  e.  , mass  is  being  pumped 
out  of  the  recirculation  region  through  the  mixing  layer  as  shown  in 
Fig.  Z lb. 


V 

j 


'Vv 


Reattachment  (Recompression)  Process 

Perfect  mass-pumping  effectiveness  is  assumed,  and  the  cutoff 
location  of  the  free  - shear- layer  region  is  consequently  ended  at  the 
rear  stagnation  point.  Furthermore,  a Nash  factor  of  unity  is  used 
which  implies 


Pp3d 

P3e 


P± 

P3 


(51) 


where  p4  is  the  pressure  that  exists  on  a wedge  of  angle  0,2  behind 
the  oblique  shock  when  a flow  with  Mach  number  M3e  exists.  Note 
that  when  M3d  > 1,  p4  = po3dy,  i.  e.  , the  value  of  po3d  is  therefore 
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The  d streamline  location, 
asymptotic  velocity  profile 

erf  t)  = 2 <t>  - l 


and  then  I}d,  the  mass-flow  integral,  can  be  determined 


Conservations  of  Mass,  Momentum 

Mass  Conservation.  Mass 
following  governing  equation 


cons 


oo  or  mc  can  be  computed  from  Eq.  50  as  shown  below 
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where  Ijd  and  Ijj  are  defined  in  Eq.  57  and  49,  respectively. 


M°glg,nju.m  Conservation,  Since  the  momentum  associated  with 
the  injection  flow  is  very  small  (less  than  4%,  as  compared  with  the 
main  flow  in  all  cold-flow  test  cases),  it  is  neglected  in  the  present 
analysis.  Therefore,  the  momentum- conservation  consideration  is 
applied  only  to  the  free  shear  layer  itself,  as  indicated  in  Eq.  49  to 
locate  the  separating  streamline. 

Energy  Conservation.  Two  forms  of  energy  transfer  are  asso- 
ciated with  the  jet  mixing  layer:  one  results  from  the  total  enthalpy 
gradient  across  the  mixing  layer,  and  the  other  results  from  mass- 
flow  entrainment.  Energy-conservation  principles  yield  the  following 
governing  equation  6 


qb  +snc  +2nd 
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(61) 


where  n£  is  the  shear  work  and  heat  transfer  across  the  mixing  layer 
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and  nd  is  the  energy  transfer  by  *he  entrained  mass  flow. 
yd 
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pu  h dv 
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The  rate  of  heat  transfer  to  the  nozzle  surface,  qb,  has  not  been  con- 
sidered in  the  present  study.  Without  knowing  the  qb  value,  the 

energy  equation  cannot  be  solved  directly,  but  it  can  be  approximated 
by  the  following  equation. 


' ob 


= T (T  + T ) 

2 Ooo  OC 


(64) 


Since  the  energy  transfer  is  not  the  important  factor  for  TVC  perform- 
ance, the  use  of  the  above  approximate  equation  should  produce  only 
a very  small  error,  if  any,  in  TVC  effectiveness  evaluation. 


Computational  Procedure 

The  analysis  is  carried  out  in  the  following  manner: 
1.  A trial  value  of  p^  is  selected. 
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FIG.  24.  Centerline  Pressure  in  CEN-17  Nozzle  at  p 
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FIG.  25.  Average  Pressure  in  CEN-17  Nozzle  at  p 
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Ignoring  the  shear  force, 
BLTVC  applications: 


Tw'  gives  the  following  equation  for 
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SIMILARITY  PARAMETER  AND  NOZZLE  DESIGN  GUIDELINES 

BLTVC  SIMILARITY  PARAMETER, 

Simulation  criteria  are  very  important  in  scaling  a BLTVC  system 
that  is  to  be  tested  in  a combined  aerodynamic -thruster  environment. 

It  is  necessary  to  maintain  geometric  scaling  and  jet-plume-to-free- 
stream  momentum  scaling,  to  match  Mach  number  and.  where  possi- 
ble, to  match  Reynolds  number  for  general  aerodynamic  similarity. 

It  is  also  necessary  to  closely  approximate  the  effective  specific-heat 
ratio  of  the  actual  main  nozzle.  Failure  to  satisfy  these  requirements 
will  preclude  matching  the  crucial  internal  nozzle  pressure  distribution 

In  AVS/dJiSirabue  t0  d“Plicate  the  energy-diffusing  parameter, 

f , c'  where  - molecular  weight,  for  the  simulation 

of  the  external  jet/free-stream  interaction. 

In  addition  to  the  considerations  noted  above,  the  results  of  the 
present  study  indicate  that  it  is  necessary  to  simulate  an  additional 
parameter,  £,  to  account  for  the  TVC  gas  interaction  with  the  nozzle 
main  flow.  C is  a two-stream  mixing -energy  parameter  that  is  a func- 
tion of  control-port  mass-flow  rate,  nozzle  mass-flow  rate,  nozzle 
total  temperature,  nozzle  molecular  weight,  missile  recovery  temper- 
ature, and  air  molecular  weight. 


Two-Stream  Mixing  Considerations 

One  method  of  extending  conventional  wind-tunnel  simulation 
requirements  to  flows  involving  two  fluid  streams  is  to  require,  in 
addition,  that  the  ratios  of  all  relevant  reference  quantities  of  the  two 
streams  be  matched,  and  that  dimensionless  mass  and  energy  param- 
eters be  included  in  the  list  of  similarity  variables.  Since  free-shear 
layers  are  nearly  always  turbulent  at  Reynolds  numbers  of  practical 
interest,  transport  by  molecular  diffusion  and  conduction  at  the  shear 
layer  may  be  neglected,  and  the  only  molecular  transport  properties 
that  must  be  included  are  those  that  influence  the  vehicle  boundary 
layer,  such  as  the  mass  and  energy  properties  of  the  control  injectant. 

-^errcPerature  Effects.  The  simulation  of  temperature  effects  can 
be  satisfied  by  matching  the  specific-heat  ratio  value  and  the  injectant' s 
mass,  momentum,  and  energy  fluxes  normalized  with  respect  to  the 
corresponding  quantity  of  the  nozzle  flow  as  shown  below. 

For  mass,  “9 

( pvA). 

°"  = w- 
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For  momentum. 


( pV  M A) . 

© = T-2- 

ipVM  A) 


For  energy 

(c  T ). 

^ = r>l  r^W- 

p o c 

where  c = (\/V-l)R  and  R is  the  gas  constant.  Subscript  j refers  to 
the  TVCPcontrol-port  injectant,  and  subscript  c refers  to  the  nozzle 
flow  at  the  separation  point. 

It  is  interesting  to  note  that  the  energy  parameter  can  be  reduced 
to  the  well-known  jet-interaction  parameter,  t.  by  the  following 
procedure: 

(c  T ) . 
p o c 

When  the  mass -flow  ratio  is  duplicated  in  the  test. 


( c„  T ) 

- P 

( c T ) 
p o c 


R T.)j 
(A  R To) 


When  V is  also  duplicated  in  the  test. 


(RTo) 

''RToK 


,.&L 
' (»  ' 


where^is  the  molecular  weig 


c 

gkt. 


T 


Therefore,  it  is  clear  that  unless  the  mass-flow  parameter,  111  , 
and  the  specific  heat  ratio,  V , are  duplicated,  matching  of  the  jet 
interaction  parameter,  t,  alone  rnay  not  entirely  satisfy  the  energy 
simulation.  Also,  the  choice  of  <?  over  t for  BLTVC  simulation  is 
obvious  due  to  its  generality.  Independent  of  the  magnitude  of  the 
temperature  effect  on  TVC  performance,  the  energy  simulation  should 
be  preserved  because  of  potential  effects  on  other  related  considera- 
tions, such  as  heat  transfer  or  motor  thrust  performance. 

Typical  £ values  for  Chandler  Evans  cold-flow  and  NWC  rocket 
static  tests  are  shown  in  Fig.  36.  Values  of  111  used  for  computing  c 
are  given  in  Fig.  37.  It  can  be  seen  in  Fig.  37  that  the  momentum 
ratio  of  the  injectant  to  the  nozzle  flow,  0,  is  very  small.  Therefore, 
the  BLTVC  does  not  appear  to  involve  a dominantly  strong  lateral-jet 
type  of  flow  interaction  under  static-jet  conditions.  , 
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Both  Fig.  36  and  Fig.  37  are  based  on  the  BLTVC  no 
rations  given  in  Fig.  38  and  Fig.  39. 


wijo  veu  Ubivo  reriormance.  j.  v u perlormance  under  repre- 
sentative hot-  and  cold-flow  conditions  for  similar  nozzle  configura- 
tions is  shown  in  Fig.  40.  The  difference  in  side  force  (or  angular 
deflection)  for  cold-flow  and  solid-propella,nt  rocket  static  tests  over  a 
range  of  chamber  pressures  appears  to  be  significant.  Therefore,  the 
hot-gas  effects  such  as  temperature,  erosion,  particle  deposition,  etc. 
on  BLTVC  performance  must  be  considered. 

As  anticipated,  the  largest  effect  appears  to  be  that  of  mass  addi- 
tion, which  serves  as  a triggering  device.  An  increase  in  the  injection 
mass  does  not  always  improve  BLTVC  effectiveness  at  a given  pQc 
value,  as  can  be  seen  in  Fig.  41.  When  the  port  size  drops  from  5/8- 
to  3 /8-inch  diameter,  the  injection  flow  area  is  reduced  by  a factor 
of  about  3,  while  the  drop  in  side  force  is  of  a lesser  order.  This 
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FIG.  38.  Ctundler-Evans  Cold-Flow  Test-Nozzle  Configuration. 
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PORT  OPENING:  ALL  THREE  HOLES 
5/8-IN.  DIAM. (PORT  SIZE) 

— • 1/2-IN.  DIAM. 

— 3/8-IN.  DIAM. 

1/4-IN  DIAM 


. i - 30 
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w 


CHAMBER  PRESSURE  (PSIG) 
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FIG.  41.  Effect  of  Mass  Addition  on  BLTVC  Effectiveness  (Cold-Flow 
Test  Results). 

NOZZLE  DESIGN  GUIDELINES 

a BI  TV?11"68.  f°r  determinin8  Port  locations  and  sizes  in  the  design  of 
fiTVC,  nozzle  afa  a oocessarj-  adjunct  to  empirical  data  on  operation 
1 fi  variotls  conditions.  For  example,  suppose  a cold-flow  thrust- 
deflection  angle  of  15  degrees  has  been  achieved  with  p = 2 500  osia 
* 5?'  a"d  l The  question  is  how  to  obtain  simlfar  performance 

Tnd\h  1§2S  lfT  thft  °Pe.ratin§  conditions  are  poc  = 1,  600  psia,  e = 50. 

, , ' ' Present  study,  a preliminary  design  guideline  was 

developed  and  is  outlined  below.  g guideline  was 

Port  Location  Determination 

interes"S‘5nCFiT°«r,VheS'  1oc/Pr  versus  <s  for  the  two  cases  of 

lines  for  two  i 1,  ® ^Ve  been  ShoWn  as  dashed  lines  and  "Olid 

lines  tor  two  different  operating  cases.  For  each  case  conditions  are 

ardetcribed^fth."  qWithT  ^ TheSe  conditi°"s  ^ obtained 
as  described  in  the  Separation  Point  section. 

ing  torr^T^8  boUndaries  such  as  A'B'C'D'  and  abed  accord- 
ing to  the  poc  limits  shown  on  the  right-hand  scale  of  Fig.  42. 

, _ th4.tP?rtS  *s  shown  in  the  lower  Portion  for  GEN-17 

h 7 .W1^^an  the  boundary  of  A'DL  Similar  TVC  performance  could 

be  achieved  for  hot-gas  conditions  if  port  locations  were  within  the 
boundary  of  ad,  as  is  also  shown  for  c = 50. 
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If  f = 70  is  required  because  the  operating  pressure  is  being  raised 
to  poc  = 2,  500  psia  and  the  boundary  is  consequently  being  changed  to 
aD,  a five-port  design  is  then  required  to  achieve  similar  performance 
because  the  operating  boundary  (ABCD)  is  larger.  Typical  port  loca- 
tions are  also  shown  for  e = 70. 

Port  Sizes 

Port  sizes  can  be  determined  according  to  the  BLTVC  similarity 
parameter,  as  discussed  in  detail  earlier  in  this  section. 

Operating  poe  Range 

It  is  desirable  to  select  a lower  poc  range  (500  to  1,  600  psia)  for 
wind-tunnel  tests  because  a higher  pressure  exists  at  the  separation 
point,  psep  (see  Fig.  43)  and  because  a shorter  nozzle  such  as  « = 50 
(see  Fig.  42)  can  be  used.  Depending  on  mission  requirements,  it 
may  be  desirable  to  select  a high  pQC  range  (500  to  2,  500  psia)  for 
flight  because  a higher  thrust  level  results,  and  longer  operating  times 
are  experienced. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


A key  objective  of  this  study  was  to  provide  analytical  tools  to 
assist  in  a largely  empirical  BLTVC  design  process.  This  has  been 
achieved  through  quantification  of  the  influence  of  dominant  BLTVC 
flow-field  phenomena.  Significant  conclusions  regarding  the  gener- 
ation of  control  forces  are  given  below. 


High  side  force  can  be  achieved  by  properly  choosing  parameters 
that  control  the  port-side  and  the  attached-side  separated  flows.  For 
the  port  side,  early  separation  is  desirable.  For  the  attached  side, 
reattachment  near  the  nozzle  aft  end  is  desired.  It  is  extremely  im- 
portant that  the  flow  reattach  on  the  attached  side  so  as  to  maintain  a 
much  lower  pressure  there  than  on  the  port  side.  Reattachment  pro- 
vides a shield  to  prevent  back-pressure  feedback,  which  results  in 
higher  pressure,  as  usually  observed  on  the  port  side.  The  cylindri- 
cal aft  body  provides  a good  mechanism  to  ensure  an  attached  flow  on 
the  attached  side.  However,  for  lower  chamber  pressures,  say 
Poc  < 1,  000  psia,  the  flow  does  not  attach,  even  with  a cylindrical 
section.  Additional  work  to  investigate  desirable  aft-section  geometry 
for  low  p0c  values  is  recommended. 

Study  of  the  separation-point  locations  and  sizes  has  provided  a 
clearer  understanding  of  how  the  TVC  gas  interacts  with  the  nozzle 
flow  and  has  helped  to  provide  a useful  guideline  for  BLTVC  nozzle 
design  and  a method  for  locating  the  TVC  ports.  The  first  port  should 
be  placed  at  or  slightly  downstream  of  the  separation  point  without  TVC. 
Failure  to  do  this  will  result  in  lower  side  force  values,  as  shown  in 
Fig.  44  and  45.  Other  effects  of  port  Location  are  shown  in  Fig.  46. 

The  effect  of  port  size  was  discussed  in  the  section  on  flow-field 
models  and  is  shown  in  Fig.  41.  The  internal  strakes  used  in  some 
designs  to  obtain  a two-dimensional-type  flow  have  pronounced  effects 
on  the  flow-field  disturbance  and  consequently  on  the  resulting  side- 
force  value.  As  can  be  seen  from  Ref.  19,  very  large,  nonuniform 
pressures  exist  in  the  vicinity  of  the  strake  leading  edge.  Thus,  small 
changes  in  location  are  important.  The  sensitivity  of  the  leading-edge 
location  may  become  even  more  critical  in  hot-gas  tests  because  of  the 
deposit  of  exhaust  products  on  the  strake  leading-edge  surfaces. 

Hence,  in  general,  the  position  of  the  strake  leading  edge  should  be 
further  downstream  for  hot-gas  tests  than  for  cold-flow  tests.  The 
shape  and  location  of  the  strake  leading -edge  effect  should  be  studied 
further,  both  analytically  and  experimentally. 
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There  is  a significant  difference  between  cold-flow  or  static 
rocket  tests  and  wind-tunnel  or  flight  tests.  As  discussed  in  detail  in 
the  section  on  flow-field  conditions,  the  ambient  flow  will  reduce  the 
flow  feedback  effect,  which  may  lower  the  pressure  on  the  port  side  — 
but  the  ambient  flow  will  also  help  flow  reattachment  on  the  attached 
side.  Consequently,  higher  side  force  may  occur  with  lower  p 
values,  as  compared  to  the  same  poc  values  in  cold-flow  tests.°C 
Therefore,  the  need  for  a cylindrical  aft  section  is  probably  less,  as 
compared  to  cold-flow  tests.  With  higher  p values,  lower  side 
forces  probably  will  be  observed  in  wind-tunnel  and  flight  tests.  This 
should  be  kept  in  mind  in  designing  a BLTVC  system  for  flight  environ- 
ments. The  combined  effects  of  ambient  flow,  altitude,  and  the  high 
temperature  of  nozzle  flow  can  best  be  evaluated  experimentally. 
Analytical  study  should  be  considered  for  optimizing  overall  perform- 
ance for  all  parameters  involved.  ^ 

The  present  study  represents  the  first  analytical  attempt  to  provide 
a method  for  BLTVC  performance  prediction,  and  cursory  assumptions 
Were  used  when  needed.  It  is  strongly  recommended  that  further 
analytical  efforts  be  conducted  to  improve  the  analytical  methodology 
and  the  prediction  of  BLTVC  performance  in  flight  environments. 


•v 


I 


■ ■ 


NWC  TP  5788 


Appendix 

AUXILIARY  CHARTS  AND  THEIR  USES 

, Subro^tines  from  an  existing  McDonnell  Douglas  computer  nroeram 

eporT^^Lr"  brte  nmVS  *"  wo*rk  covered  b ■ 

Y P°  ‘ u , to  obtain  "Pmencal  flow-field  solutions.  Specifically 

assembled  to  predict  local  flow-field  solutions' 
b arg.®  variation  of  flow  parameters  such  as  M T n v 
Tooo.\o..  and  port  size.  Thirty-six  charts  (FigC’  47  through  82) 
were  constructed  to  treat  the  following  parametric  variations?  ’ 

X K / 


M 


(at  separation  point)  = 3.  5,  4,  and  4.  5 


Toc  = 3,000,  4,000,  and  5,  500#R 
Ye  = 1.  2 


=1.2  and  2 

To«=  5°o°R 


= 1.4 


3,  4,  6,  and  8 


Y« 

P«>/Pc 
m/pcuc  = 0 to  0.  2 

Figures  47  through  55  present  pb/p  as  a function  of  m/pm. 
e Ph'P«  ls  an  implicit  function  of  m/n  n . 


Sin  . through  55  present  pb/p  as  a function  of  m/pm 

obtafn^XLm  lmpl‘Ci‘  fUnC“°n  “—.ion  is  neeSeS'to 

UnedAaiTono1ws?mPUtati°nal  procedure  that  asas  ‘hese  charts  is  out- 

detem  ouL^iATAe^aYthe  ^'eparaflor/poim!^*311^6  va‘ues’ 

o/u/Aa)PCfr„mC,step  1?°  ^ P"/Pc  ValUeS  baSed  °"  the  value 

until goF/ig’  47  thf°U®h  55  and  83  through  86,  iterate  m/pcuc 
until  pb/p^  is  satisfied.  Note:  rh  = (pv)w.  Pcc 

Use  pb/Pee)  to  construct  a pressure-distribution  curve  for  the 
port  side. 

Nomenclature  is  given  in  Fig.  83.  Several  auxiliary  charts  for 
expediting  computation  are  given  in  Fig.  84,  85,  and  86! 
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FIG.  85.  Nozzle  Mass  Flow  Rate  Versus  Chamber  Pressure, 
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I Naval  Ordnance  Station,  Indian  Head  (Code  FS,  A.  T.  Camp) 

1 Naval  Undersea  Center,  San  Diego  (Code  133) 

3 Naval  Surface  Weapons  Center,  White  Oak 

Code  33.  H.  Heller  (I ) 

Code  C.  C.  W.  Bernard  (1) 

Code  GW  (I ) 

1 Naval  Intelligence  Support  Center  Liaison  Officer  (LNN) 

4 Army  Armament  Command  ( AMSAR-SF) 

4 Picatinny  Arsenal  (AMD.  Concepts  Branch) 

4 Air  Force  Systems  Command,  Andrews  Air  Force  Base 

PL(1) 

DLPI(l) 

DLW(I) 

SDW(I) 


» 


2 Air  Force  Aero-Propulsion  Laboratory.  Wright-Patterson  Air  Force  Base 
R.IA  ( I ) 

STINSO  ( I ) 

5 Air  Force  Armament  Laboratory.  Eglin  Air  Force  Base 
DLD ( I ( 

DLI(I) 

DI.0  ( I > 

DLQ(i) 

DIAL  C APT  Aden  ( I ) 

I Air  Force  Rocket  Propulsion  Laboratory,  Edwards  Air  Force  Base  (Code  MKP) 

I Foreign  Technology  Division,  Wright-Patterson  Air  Force  Base  (Code  PDXA.  James  Woodard) 
I Wright-Patterson  Air  Force  Base  (Code  XRHP) 

1 Defense  Advanced  Research  Projects  Agency,  Arlington 

2 Defense  Documentation  Center 

I Purdue  University.  School  of  Mechanical  Engineering.  Lafayette.  It.d.  (Cecil  F.  Warner) 
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